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Integration  of  an  Individual-Based  Fish 
Bioenergetics  Model  into  a  Spatially  Explicit 
Water  Quality  Model  (CE-QUAL-ICM) 

by  P.  Soupy  Dalyander  and  Carl  F.  Cerco 


PURPOSE:  This  technical  note  presents  the  results  of  incorporating  a  fish  bioenergetics  module 
into  CE-QUAL-ICM,  a  spatially  explicit  eutrophication  model.  In  addition  to  fish  consumption 
of  algae,  zooplankton,  and  detritus,  fish  biomass  accumulation  and  recycling  to  the  water  column 
are  explicitly  accounted  for.  Schools  of  fish  are  tracked  individually,  allowing  for  spatial  resolu¬ 
tion  of  their  effects  on  phytoplankton  and  nutrient  loading.  These  developments  allow  for  invest- 
tigations  of  the  impact  of  planktivorous  fish  on  water  quality  through  consumption  of  algae  and 
nutrient  export  via  mass  sequestering  during  growth.  The  formulation  details  of  the  new  model 
system  are  outlined  in  this  technical  note. 

BACKGROUND:  As  early  as  the  1960’s-70’s  intense  eutrophication  events  have  occurred  in 
sensitive  estuaries,  with  increasing  anthropogenic  sources  of  nitrogen  (N)  and  phosphorus  (P) 
leading  to  excessive  algal  growth,  degradation  of  water  quality,  and  a  general  decline  in  eco¬ 
logical  health.  Most  prior  attempts  to  address  estuarine  eutrophication  have  focused  on  reducing 
nutrient  loads;  however,  cost  and  difficulty  in  implementation  have  limited  the  effectiveness  of 
such  an  approach.  An  alternate  approach  under  consideration  is  the  revitalization  of  primary  con¬ 
sumers,  including  planktivorous  fish,  which  can  potentially  reduce  the  quantity  of  algae  and 
remove  nutrients  from  estuarine  systems  via  sequestering  during  growth.  The  model  formulation 
presented  in  this  technical  note  is  designed  to  provide  a  way  to  quantify  the  potential  impact 
managing  planktivorous  fish  may  have  on  the  water  quality  of  an  ecosystem.  A  fish  bioener¬ 
getics  module  based  on  the  Wisconsin  fish  model  (Hanson  et  al.  1997)  was  coupled  into  the 
spatially  and  temporally  explicit  water  quality  model,  CE-QUAL-ICM.  Fish  are  modeled  as  indi¬ 
vidual  schools,  allowing  spatial  and  temporal  resolution  of  water  quality  effects,  and  nutrient 
feedback  through  egestion,  excretion,  and  mortality  is  explicitly  modeled. 

EUTROPHICATION  MODEL:  CE-QUAL-ICM  is  a  three-dimensional  eutrophication  model 
that  operates  by  enforcing  conservation  of  mass  on  a  set  of  state  variables  in  each  cell  of  a  model 
grid.  The  state  variables  each  belong  to  a  set  of  categories  that  may  be  described  as  physical 
variables  (salinity,  temperature,  etc.)  and  those  in  elemental  cycles  including  carbon,  nitrogen, 
phosphorus,  silica,  and  oxygen  (Cerco  and  Noel  2005;  Cerco  and  Cole  1993).  In  the  current  for¬ 
mulation,  active  constituents  are  temperature,  salinity,  fixed  inorganic  suspended  solids,  three 
algal  groups  (nominally  representing  blue-green  algae,  spring  diatoms,  and  green  algae),  two 
zooplankton  groups  (representing  micro-  and  mesozooplankton),  ammonium,  oxidized  nitrogen 
(nitrite/nitrate),  dissolved  oxygen,  and  various  labile  and  refractory  forms  of  dissolved  and  parti¬ 
culate  carbon,  nitrogen,  phosphorus,  and  silica.  Sources  to  the  model  are  meteorological  inputs 
and  loading  sources  simulating  riverine  input;  these  sources  are  time-variant  to  capture  the 
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seasonal  shifts  in  nutrient  loading.  The  full  details  of  the  CE-QUAL-ICM  model  formulation 
may  be  found  elsewhere  (Cerco  and  Cole  1993;  Cerco  and  Noel  2005). 

FISH  POPULATION  DYNAMICS:  Within  CE-QUAL-ICM,  fish  are  modeled  as  discrete 
“super-individuals,”  representative  of  groups  of  individuals  moving  throughout  the  model  grid. 
Parameters  governing  fish  physiology  are  prescribed  by  “class,”  which  may  represent  a  specific 
age,  size,  or  other  sub-group  of  fish.  Hereafter  super-individuals  will  be  referred  to  as  “schools,” 
noting  that  the  difference  between  a  super-individual  and  an  actual  school  is  the  more  diverse 
nature  of  the  latter.  A  school  may  either  be  totally  uniform  in  nature,  consisting  of  a  single  class 
of  fish,  or  it  may  consist  of  fractional  components  of  multiple  classes,  each  represented  by  an 
individual  whose  physical  characteristics  (length,  weight,  composition,  etc.)  are  identical  to  all 
those  members  within  the  school  of  the  same  class. 

Fish  movement  is  governed  by  an  external  control  file  which  prescribes  the  initial  number  of  fish 
in  a  school,  the  fractional  division  between  classes,  and  the  cell  location  at  any  given  point  in 
time.  This  file  is  generated  externally,  and  may  either  be  done  manually  or  through  the  use  of 
algorithms  to  approximate  fish  migration  with  time. 

Fish  mortality  within  a  school  is  accounted  for  within  the  model  (variables  used  in  all  equations 
are  listed  in  Appendix  A): 

AA  =  ANstv  +  A Nmf  +  A Nprd  +  A Nfsh  (1) 

The  total  number  of  individuals  lost  (AN)  within  a  time-step  (At)  is  the  sum  of  the  number  lost 
due  to  starvation  (ANstv)  and  suffocation  (ANsuf),  accounted  for  explicitly  as  functions  of  fish 
condition  and  model  environment  at  any  given  time-step,  along  with  the  number  lost  to  other 
natural  causes  (ANprd,  typically  predation)  and  those  caught  by  the  fishery  (ANfSh).  The  number 
lost  due  to  any  specific  process  is  a  function  of  the  mortality  rate  for  that  process,  e.g.  for 
starvation: 


^Nstv  ~  N,  ■  Mstv  ■  At  (2) 

M  is  mortality  (fish/s)  and  Nt  is  the  number  of  individuals  (fish)  within  a  given  school.  Starvation 
and  suffocation  were  accounted  for  based  on  models  developed  for  Chesapeake  Bay  anchovy 
(Adamack  2007).  If  the  wet  weight  of  the  representative  individual  for  a  school  drops  below 
70  percent  of  the  ‘healthy’  weight  for  that  individual’s  length,  a  daily  mortality  rate  of 
0.1  fish/day  is  introduced  (converted  to  instantaneous  mortality  rate  in  fish/s  as 
Mstv  =  DFST /  (24-60-60),  where  DFST  is  the  daily  mortality  rate).  ‘Healthy’  wet  weight  (WW,  in 
g)  is  calculated  from  empirically  determined  length  (L,  in  mm)  to  weight  relationships,  specific 
to  fish  species,  and  is  presently  fonnulated  as: 

WWindMalthy  =  AL-LBL  (3) 

where  AL  and  BL  are  empirically  determined  parameters.  Suffocation  mortality  is  calculated  as 
(Adamack  2007): 
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DFSU  =  0.093487  +  70. 1 1 894  •  (in  [DOcel,]f  if  [ DOcell ]  <  1  mg // 

=  0  if[DOcell]>lmg/l 

[DOceii]  is  the  concentration  of  dissolved  oxygen  {mg/I).  Predation  and  fishing  mortality  are 
taken  from  values  reported  in  the  literature. 

FISH  BIOENERGETICS  MODEL:  Bioenergetics  models  rely  on  conservation  principles,  with 
the  amount  of  energy  a  fish  can  devote  to  growth  being  the  difference  between  consumption  in 
the  form  of  food  and  the  sum  of  life  process  expenditures,  including  respiration,  specific  dyna¬ 
mic  action,  egestion,  and  excretion.  The  rate  of  wet  weight  growth  can  be  calculated  as: 

^-{C’-KR'+S'MF'+U’)]}-^-  (5) 

01  npred 

dWW/dt  is  the  rate  of  wet  weight  growth  (g/s);  Epred  is  the  energy  density  of  the  fish  (J/g).  The 
use  of  the  prime  symbol  denotes  these  are  rates  of  energy  uptake  of  consumption  (C'),  respiration 
(R'),  specific  dynamic  action  (S'),  egestion  (F'),  and  excretion  (U'),  as  opposed  to  a  total  quantity. 
Fish  energy  density  is  calculated  as  a  linear  function  of  representative  individual  weight: 

Epred,t  =  EPRD  ■  WWind  t  (6) 


EPRD  is  an  empirically  detennined  predator  energy  density,  variant  by  species.  Growth  in 
weight  can  subsequently  be  converted  to  growth  in  length  via  the  aforementioned  empirical 
length/weight  relationship,  noting  that  fish  cannot  shrink.  At  each  time-step,  the  length  of  a 
school’s  representative  individual  is  calculated  from  its  new  weight  and  compared  to  its  length  at 
the  previous  time-step;  if  the  new  length  is  greater  than  the  old,  the  fish  grows  to  this  value,  but 
if  not,  the  fish  retains  its  prior  length  and  begins  to  have  a  weight  deficit  that  could  potentially 
lead  to  starvation  mortality  as  previously  described. 

Consumption  is  based  on  a  foraging  model  modified  from  one  developed  for  Chesapeake  Bay 
menhaden  (Luo  et  al.  2001).  Presuming  the  number  of  fish  within  a  cell  relative  to  the  cell 
volume  and  length  of  time-step  is  not  a  limiting  factor,  each  individual  within  a  school  has  a 
volumetric  clearing  rate: 

v'  =  f(DO)  ■  gap{L)  -u(T,L)-  eff  (L)  (7) 

V's  is  the  volumetric  clearing  rate  (m3/s)  and  is  based  on  the  mouth  gap  area  (gap,  in  nr),  swim¬ 
ming  speed  (u,  in  m/s),  an  empirically  derived  function  of  dissolved  oxygen  concentration  within 
the  cell  the  school  occupies  (f(DO)),  and  filtration  efficiency  (eff(L),  unitless).  The  dependence 
on  dissolved  oxygen  is  based  on  physiological  fish  response  to  stress  (Luo  et  al.  2001): 


f{DO) 


_ 1 _ 

1  +  exp  (-2. 1972  •  [DO]  +  6.5916) 


(8) 
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[DO]  is  the  concentration  of  dissolved  oxygen  (mg/1).  The  mouth  gap  area  is  a  function  of  body 
length  (Luo  et  al.  2001): 

gapi^L)  =  0.2586e-5  -(L)1  79767  (9) 


Swimming  speed  in  some  filter  feeding  species  is  a  function  of  body  length  and  food  concentra¬ 
tion,  with  fish  exhibiting  one  characteristic  swimming  rate  (in  body  lengths/s)  when  chlorophyll 
concentration  exceeds  a  threshold  to  trigger  ‘feeding’  behavior  and  a  second,  slower  rate  when 
chlorophyll  does  not  exceed  this  threshold  (Durbin  and  Durbin  1975).  In  the  present  fonnulation, 
based  on  observations  of  Atlantic  menhaden,  swimming  velocity  is  calculated  as: 

u  =  VELF  ■  if  Conrm  >  FTHR 
1000  c 

=  VELNF  ■  if  Con,  .,.  <  FTHR 

1000  •  (10) 

VELF  and  VELNF  are  the  swimming  velocities  when  feeding  and  non-feeding,  respectively 
(fish  body  length/s),  ConCHL  is  the  concentration  of  chlorophyll  (g/m  ),  and  FTHR  is  the  feeding 
threshold  (in  g/m  ).  Modification  would  be  necessary  when  applied  to  a  species  with  a  different 
observed  behavior.  Swimming  velocity  is  consistent  with  the  migration  velocity  used  in  tracking 
school  location  in  that  a  school  cannot  migrate  an  unreasonably  large  distance  compared  to 
typical  menhaden  swimming  velocities,  but  the  two  are  not  directly  coupled.  A  direct  link  is  not 
required  since  migration  velocity  captures  directionally  specified  movement  on  relatively  large 
temporal  and  spatial  scales  (hours  and  km),  whereas  moment-to-moment  fish  swim  velocity  and 
direction  on  the  scale  of  the  water  quality  time-step  (minutes)  may  be  much  more  sporadic. 

Filtration  retention  efficiency  may  also  be  a  function  of  fish  length.  In  the  current  formulation, 
the  relationship  is  a  sigmoid  response  curve  developed  for  Atlantic  menhaden  (Luo  et  al.  2001): 

eff(L)  =  t - 7 - — - rv  (11) 

v  7  (l  +  exp (-0.052781  \-L  +  2.96973)) 


In  formulations  with  application  to  other  species  of  fish,  modification  of  this  empirically  deter¬ 
mined  relationship  may  be  required. 

Within  the  model  cells,  it  is  possible  that  the  volumetric  clearing  rate  by  schools  will  need  to  be 
limited  so  that  the  volume  filtered  within  the  time-step  does  not  exceed  the  cell  volume.  In  this 
case,  the  fraction  of  potential  volumetric  clearing  rate  by  an  individual  is  compared  to  the  total 
volumetric  clearing  rate  of  all  individuals  within  the  cell.  This  value  is  multiplied  by  the  volu¬ 
metric  clearing  rate  to  turn  over  the  entire  cell  during  a  time-step  to  determine  the  maximum 
clearing  rate  for  that  individual.  Fish  mortality  is  accounted  for  prior  to  feeding,  therefore  the 
number  of  individuals  in  a  school  is  the  value  at  the  end  of  the  time-step: 


V' 

s,  max 


AV 

At 


K 


(12) 
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In  this  equation  AV  is  the  total  cell  volume  (m3);  At  is  the  length  of  the  time-step  (s);  V's  is  the 
volumetric  clearing  rate  of  an  individual  within  the  school  as  defined  previously;  and  NSjt+At  is 
the  number  of  individuals  at  that  end  of  the  time-step  (e.g.,  after  mortality).  If  the  combination  of 
fish  population,  cell  size,  and  model  time-step  are  such  that  this  alternative  calculation  is  widely 
employed,  growth  or  oxygen  intake  will  become  unrealistically  restricted,  and  the  grid  cell  size 
or  the  time-step  will  need  to  be  modified. 

Energy  consumption  is  a  function  of  volumetric  clearing  rate,  filtration  efficiency,  and  the 
energy  concentration  within  the  grid  cell,  calculated  from  the  mass  density  of  food  sources  and 
their  energy  density.  Model  fish  consume  food  at  all  times,  causing  unrealistically  large  overall 
growth  compared  to  biological  fish,  which  only  consume  food  during  a  portion  of  their  waking 
hours.  Therefore,  a  feeding  fraction  is  introduced,  which  limits  food  consumption  at  all  times  to  a 
fraction  of  possible  total  consumption.  This  feeding  fraction  is  selected  empirically  to  calibrate 
fish  growth  patterns  to  those  observed  in  the  field. 

The  three  types  of  prey  considered  within  the  model  are  phytoplankton,  mesozooplankton,  and 
detritus.  The  total  algal  concentration  is  the  sum  of  the  three  model  algae  species  (Conphyt0,  in 
gC/nr3);  retention  efficiency  is  independent  of  prey  source  in  the  model. 

Detrital  food  sources  are  also  considered.  Detritus  is  handled  by  the  model  as  individual  ele¬ 
mental  forms  of  particulate  matter,  the  sum  of  which  is  taken  as  the  dry  weight  of  detritus: 

Condet  =  LPOC  +  RPOC  +  LPON  +  RPON  +  LPOP  +  RPOP  +  PS  (13) 

-3 

Condet  is  the  total  (dry)  concentration  of  particulate  constituents  (g/m  ),  including  labile  and 
refractory  particulate  organic  carbon  (LPOC  and  RPOC,  in  gC/m  ),  labile  and  refractory  particu¬ 
late  organic  nitrogen  (LPON  and  RPON,  in  gN/  m3),  labile  and  refractory  particulate  organic 
phosphorus  (LPOP  and  RPOP,  in  gP /  m  ),  and  particulate  silica  (PS,  in  gSi/  m  ).  Concentration 
of  wet  weight  of  detritus  is  calculated  assuming  a  dry  weight  to  wet  weight  ratio  of  0.2  (Rippetoe 
1993),  and  a  fixed  energy  density  set  in  the  control  file  is  used  to  convert  this  wet  weight  concen¬ 
tration  to  energy  concentration.  These  calculations  yield  a  very  approximate  value  of  energy  con¬ 
sumption  from  detrital  sources  for  menhaden,  given  that  actual  detritus  energy  density  has  a 
complex  dependence  on  composition  that  could  not  be  captured  using  the  model  constituents 
(Tenore  1981). 

Combining  phytoplankton,  zooplankton,  and  detrital  food  sources,  the  total  rate  of  energy  con¬ 
sumption  becomes: 


C'  =  FFRC  ■ 


Con 


phyto 


CtWP 


EPLK  +  - 


LZ 

CtWZ 


■EZOO  +  - 


Con 


det 


dwwwa 


■EDET 


det 


■V' 


(14) 


LLRC  is  the  aforementioned  empirically  detennined  feeding  fraction  (unitless).  Algal  model 
units  are  grams  of  carbon/m  ,  but  energy  density  values  (EPLK)  are  in  joules  per  gram  of  wet 
weight,  therefore  a  conversion  factor  of  carbon  to  wet  weight  ratio  (CtWP,  in  gC/gWW)  is 
required,  and  is  taken  as  0.1  (Peters  and  Downing  1984).  Similarly  for  zooplankton, 
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CtWZ  =  0.04344,  based  on  a  gDW/gWW  ratio  of  0.1086  (Durbin  and  Durbin  1998)  and  a 
gC/gDW  ratio  of  0.4. 

Energy  is  lost  due  to  physiological  processes  including  respiration,  specific  dynamic  action, 
excretion,  and  egestion,  which  are  modeled  using  a  methodology  following  the  Wisconsin  Fish 
Model  (Hanson  et  al.  1997).  Rate  of  energy  loss  to  respiration  is  detennined  from  mass  con¬ 
sumption  of  oxygen.  R'o  is  calculated  from  physiological  parameters  specific  for  a  given  fish 
species,  a  temperature  dependence  function,  and  an  activity  multiplier. 


R'  =  R'o-Eo 

(15) 

r’0  =  resa-ww,MS£  ■f(r)-JCT-(m%6400 ) 

(16) 

R’o  is  the  rate  of  oxygen  consumption  (g/s);  Eo  is  an  oxycalorific  coefficient  (J/g  O2);  RESA  and 
RESB  are  the  fish-specific  physiological  parameters  on  an  annual  basis;  f(T)  is  a  temperature 


dependency  function,  and  ACT  is  an  activity  multiplier.  The  temperature  dependence  of  oxygen 
consumption  is  calculated  as: 

f(T)  =  vx  ■exV(x\\-v)) 

(17) 

V  =  (RTM  -T)/ (RTM  -RTO) 

(18) 

X  =  [z2  •  ( 1  +  (1  +  40/7)°'5  J  j ^400 

(19) 

Z  =  In  (RQ)  •  ( RTM  -  RTO ) 

(20) 

Y  =  In  (RQ)  ■  (RTM  -  RTO  +  2) 

(21) 

T  is  temperature  (°C)  and  RQ,  RTM,  and  RTO  are  fish-specific  physiological  parameters.  Since 
swimming  velocity  depends  on  the  chlorophyll-triggered  feeding  state,  a  similar  behavior  depen¬ 
dence  is  introduced  into  both  the  activity  multiplier  and  the  oxycalorific  coefficient  based  on 
empirical  observation,  with  E0,  feed  =  13388.8  J/g  O2;  Eo,  nonfeed  =  13723.5  J/g  O2  (Durbin 
and  Durbin  1983).  Activity  multipliers  during  both  feeding  and  non-feeding  states  are  set  in  the 
fish  control  file. 

The  bioenergetics  model  assumes  unlimited  oxygen  availability,  which  may  not  be  true.  If 
oxygen  is  limiting  (i.e.,  if  the  fish  would  consume  more  oxygen  in  a  time-step  than  is  present  in 
the  grid  cell),  the  volumetric  clearing  rate  previously  calculated  is  used  to  determine  the  actual 
rate  of  oxygen  mass  consumption: 

*o-K-.\DO]  02) 
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Egestion  energy  loss  rate  is  calculated  as  a  percentage  of  consumption,  and  excretion  and  speci¬ 
fic  dynamic  action  losses  are  calculated  as  a  percentage  of  rate  of  energy  assimilated  (e.g.,  con- 


sumption  minus  egestion). 

S'  =  SDA\C'-F') 

(23) 

F'  =  FA-C' 

(24) 

U'  =  UA\C'-F ') 

(25) 

SDA,  FA,  and  UA  are  empirically  determined,  species-specific  factors. 

Fish  Composition:  Fish  bioenergetics  calculations  are  based  on  energy  conservation,  but 
within  the  model  environment,  constituent  element  mass  (e.g.,  carbon,  nitrogen,  phosphorus, 
silica)  must  also  be  conserved.  In  the  model,  an  empirical  value  of  fish  composition  for  each  of 
the  constituent  elements  was  specified  for  each  class,  fixing  the  ideal  fraction  of  each  element  as 
a  function  of  total  wet  weight.  Under  certain  environmental  conditions,  it  is  possible  that  con¬ 
sumption  of  an  individual  element  relative  to  growth  detennined  by  energy  consumption  is  insuf¬ 
ficient  to  maintain  this  target  fraction,  therefore  elasticity  is  allowed  wherein  a  fish  can  incur  an 
elemental  deficit,  and  subsequently  retain  higher  fractions  of  that  element  when  it  is  in  abun¬ 
dance  to  regain  the  target  composition.  The  amount  of  a  given  element  a  fish  needs  to  retain  at 
each  time-step  is  calculated  as: 

Need E  =  WWl+At  ■  D  WWW  ■  FED  W  -  WWr  ■  D  WWW  ■  aFED  Wt  (26) 


This  amount  (Needs)  is  a  function  of  the  dry  weight  to  wet  weight  ratio  (DWWW,  unitless),  the 
wet  weight  before  and  after  the  time-step  (WW,  g),  the  target  fraction  of  that  element  (FEDW, 
unitless),  and  the  actual  fraction  of  that  element  at  the  previous  time-step  (aFEDWt,  unitless).  It 
should  be  emphasized  that  the  elasticity  in  the  model  does  not  mimic  realistic  changes  in  fish 
composition;  however,  it  allows  for  conservation  of  mass  within  the  model  and  prevents 
unrealistic  individual  element  consumption  requirements  from  being  imposed.  Elemental  compo¬ 
sition  could  theoretically  drop  to  unrealistically  low  levels,  but  such  a  deficit  would  indicate  a 
problem  with  the  parameterization  of  the  model  environment  or  the  fish,  as  a  real  fish  is  still 
subject  to  elemental  mass  availability  and  conservation. 

Coupling  to  Water  Quality  Model:  To  conserve  mass  on  the  model  grid,  prey  constituent 
concentrations  must  decrease  with  consumption,  and  fish  outputs  due  to  respiration,  egestion, 
and  excretion  must  return  to  the  water  column.  In  addition,  mass  from  fish  that  die  of  natural 
mortality  should  return  in  the  grid. 

For  each  algal  group,  a  school’s  consumption  rate  is  calculated  from  the  volumetric  clearing  rate, 
the  number  of  individuals  in  the  school,  and  the  original  concentration: 


dt 


-V;-Nl+At-Bk 


(27) 
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Bk  gives  the  concentration  of  algae  group  k  (gC/nr).  Similarly  for  mesozooplankton: 

8J^-  =  -V;-N^-LZ  (28) 

at 

LZ  gives  the  concentration  of  mesozooplankton  (gC/nr).  Conservation  of  other  constituents  is 
relatively  more  complicated,  in  that  fish  can  either  be  solely  producers  (e.g.,  ammonium)  or  both 
consumers  and  producers  (e.g.,  particulate  matter).  For  each  element  considered  within  the 
model  (carbon,  nitrogen,  phosphorus,  silica),  a  recycling  rate  is  calculated  as  the  total  rate  at 
which  the  element  is  returned  to  the  water  column  by  the  fish. 

Carbon  mass  balance  dictates  the  growth  rate  in  carbon: 

G<c=C'c-(R'c+F^U<c)  (29) 


G'c  is  the  growth  rate  in  carbon  (gC/s)  and  is  a  function  of  carbon  uptake  (C  c,  in  gC/s),  respire- 
tion  ( R'c ,  in  gC/s),  egestion  (F'c,  in  gC/s),  and  excretion  rates  (U'c,  in  gC/s).  Carbon  uptake  is 
governed  by  fish  consumption,  a  function  of  the  actual  volumetric  filtering  rate  of  the  school  and 
the  original  concentration  of  the  carbon  containing  prey  sources  (algae,  mesozooplankton,  labile 
and  refractory  organic  carbon): 


r  3 


C  =  V  ■  N, 


Y  Bk  +  LZ  +  LPOC  +  RPOC 


\k= 1 


(30) 


Inorganic  carbon  is  not  explicitly  accounted  for  in  the  model;  however,  the  loss  of  carbon  due  to 
respiration  must  be  considered.  The  rate  of  carbon  emission  due  to  respiration  is  calculated  as: 

R'c  =  R'0/AOCR  (31) 


R'c  is  the  rate  of  carbon  emission  (gC/s)  and  is  a  function  of  the  rate  of  oxygen  consumption  and 
the  oxygen-to-carbon  ratio  in  respiration  (AOCR,  in  gCF/gC).  Using  the  same  simple  model  for 
aerobic  respiration  as  in  the  water  quality  model: 

C6Hu06  +  6  02  6  C02  +  6  H20  (32) 


Respiration  of  6  moles  of  dissolved  oxygen  produces  6  moles  of  carbon  dioxide,  which  converts 
to  an  AOCR  of  2.67  g  02  g'!C. 

The  ideal  carbon  recycling  rate  (CRRATE,  in  gC/s/fish)  is  therefore: 


( 


CRRATE  =  Fc+U'c=C'c-G'c -R'c  =  C'c-Needc  - 


R' 


AOCR 


(33) 


If  the  calculated  recycling  rate  were  to  fall  below  zero,  the  fish  would  (unrealistically)  absorb 
nutrients  directly  from  the  water  column;  therefore,  the  actual  carbon  recycling  rate  is  calculated 
as: 
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CRRATE  =  max 


Cc  -  Needc 


'  Rq  " 

kAOCR , 


The  total  change  of  any  non-particulate  constituent  j  is  calculated  as: 

dCj  Nr+At  ■  CRRA  TE  ■  FCF  +  Nt  ■  Mnat  ■ WWt-D  WWW  ■  FCD  Wt  ■  FC  DM 

~7T~  av 


(34) 


(35) 


The  distribution  of  recycled  carbon  to  a  model  constituent  (denoted  as  ‘j’;  e.g.,  refractory  parti¬ 
culate  organic  carbon,  labile  dissolved  organic  carbon,  etc.)  by  a  living  fish  (FCjF)  is  governed 
by  fractions  fixed  as  input  parameters  to  the  model.  Note  that  fish  consume/produce  discrete 
masses  of  constituents,  but  constituent  quantities  are  accounted  for  as  volumetric  concentrations. 
Additional  carbon  recycling  occurs  due  to  natural  mortality  (starvation,  suffocation,  and  preda¬ 
tion,  the  sum  of  which  results  in  a  total  natural  mortality  rate,  Mnat);  in  contrast,  fishing  mortality 
is  a  sink  from  the  system.  The  distribution  of  recycled  carbon  from  mortality  is  controlled  by 
fractions  input  to  the  model  (FCjFM). 

For  particulate  matter,  losses  to  fish  consumption  are  also  included: 


SC,  Nnu-  CRRA  TEf  ■  FC,f  +  N,-Mmt-WWrD  WWW  ■  FCDW,  V's  ■  [y] 

dt  AV 


(36) 


3 

In  this  equation,  [j]  is  the  concentration  of  constituent)  (in  g/m ).  Note  that  fish  die  before  they 
eat,  therefore  the  “living”  fish  recycling  rate  is  multiplied  by  the  new  population,  whereas  the 
recycling  due  to  mortality  is  multiplied  by  the  old  population  and  mortality  rate.  The  recycling 
fractions  for  a  specific  fish  may  be  estimated,  for  example,  from  empirically  determined  values 
for  that  species  or  a  similar  one.  The  new  fish  carbon  fraction  is  calculated  from  the  actual  con¬ 
sumption  and  recycling  rates  as: 


FDCW, 


(\VWt  ■  DWWW ■  FDCW,  +  A t\Cc-Rc- CRRATE f )) 
DWWW-Wt+At 


(37) 


For  the  nitrogen  constituents,  respiration  is  not  a  factor,  therefore  nitrogen  growth  rate  is  gov¬ 
erned  by: 

G'„=Ck-(FI+U'„)  (38) 


Within  the  eutrophication  model,  algal/zooplankton  nitrogen  is  accounted  for  as  a  fraction  of 
each  species’  carbon  concentration.  Nitrogen  uptake  is  governed  by  this  fixed  ratio  and  the  nitro¬ 
gen  in  the  particulate  matter  consumed: 


C'  =  V'  ■ 

V  s 


IK  •  Anck  )  +  LZ ■  ANClz  +  [LPON]  +  [RPON] 


k= 1 


(39) 
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Anck  is  the  algal  nitrogen  to  carbon  ratio  (gN/gC);  ANClz  is  the  zooplankton  nitrogen  to  carbon 
ratio  (gN/gC).  Following  the  calculations  for  carbon,  the  required  nitrogen  intake  to  maintain/ 
achieve  the  target  composition  is  calculated  as: 

NeedN  =  WWt+At  ■ DWWW  ■  TFNDW  -  WWt  ■  DWWW  ■  FNDWt  (40) 

Nitrogen  recycling  to  the  water  column  is  therefore: 

NRRATEf  =  F'n+U'n=C'n-G'n  =  max [C'N  - NeedN, 0]  (41) 


And  the  rate  of  recycling  of  dissolved  nitrogen  becomes: 

SNj  _  Nt+ A,  ■  NRRA  TEf  ■  FNjf  ~  N,  ■  Mna,  '  W,  ■  D  WWW  ■  FND  W,  ■  FNjfM 

dt  ~  AV 

Whereas  for  particulate  constituents  the  rate  of  change  becomes: 

8NJ  AU- NRRATEf  -FN/  +  N,  -Mnal  ■ WW ] -DWWW-FNDW, -FN ,/M - Nni,  ■¥;■[,] 
dt  AV 


(42) 


(43) 


The  recycling  fraction  of  nitrogen  to  each  model  constituent  from  egestion/excretion  is  again 
based  on  empirical  observation.  The  new  fish  composition  is  calculated  from  consumption  and 
recycling  rates  as: 


FNDWI+At 


(W,  ■ DWWW  ■  FNDW,  +  A t-(CN-Rc-  NRRATEf )) 
DWWW -Wn  Al 


(44) 


Calculations  for  phosphorus  constituents  are  identical  to  those  for  nitrogen  constituents. 

The  fraction  of  silica  within  fish  is  taken  as  zero,  therefore  all  silica  is  recycled  immediately  to 
the  two  model  constituents  in  a  50/50  split: 


GDS  FDSf 


dt 

AV 

dPS  _ 

f  FDSf 

dt 

{  AV 

1  v  l+Al  '  s 


£(^<vs,)+[re] 


k= 1 


■N, 


t+At 


Z(^(.S,)  +  [PS] 


k=l 


-N^-K-lPS] 


(45) 

(46) 


[DS]  and  [PS]  are  the  concentrations  of  dissolved  and  particulate  silica,  and  the  other  variables 
within  the  equation  follow  the  convention  previously  described  for  carbon,  nitrogen,  and  phos¬ 
phorus.  Rate  of  change  of  concentration  of  dissolved  oxygen  is  related  to  the  mass  of  oxygen 
consumed  per  fish,  the  population,  and  the  cell  volume: 
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dDO 

dt 


-N, 


t+At 


AV 
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(47) 


Several  assumptions  and  simplifications  are  inherent  to  this  estimate  of  nutrient  recycling  from 
the  fish  to  the  model  environment.  Both  mortality  and  the  life  processes  of  respiration,  egestion, 
and  excretion  are  treated  as  instantaneous  processes,  with  mass  returning  to  the  water  column  in 
the  same  time-step  as  death/food  consumption.  In  addition,  elemental  mass  state  must  be  simpli¬ 
fied  to  fit  within  the  confines  of  the  water  quality  model  constituents. 

MODEL  IMPLEMENTATION:  The  fish  bioenergetics  module  is  turned  on  and  off  using  a 
control  card  added  to  the  water  quality  model  control  file,  at  present  hard-wired  to  be  named 
“wqm  con.npt”  (Appendix  B;  note  that  only  a  portion  of  the  water  quality  model  is  shown  to 
illustrate  the  proper  placement  of  the  fish  card). 

A  separate  control  file,  “fbm  con.npt,”  governs  the  implementation  of  the  fish  bioenergetics 
module  (Appendix  C).  The  first  card  is  the  TITLE  card,  and  consists  of  six  lines  of  text  used  to 
make  notes  regarding  the  model  run  (these  lines  do  not  impact  actual  code  implementation). 

The  next  control  card,  “JDO,”  indicates  the  julian  day  offset  of  the  fish  migration  in  relation  to 
the  julian  day  within  the  model  run.  If  the  fish  migration  file  is  referenced  to  the  same  start  time 
as  the  model  run,  this  card  should  be  set  to  zero.  If  the  fish  migration  is  referenced  to  some  other 
point  (for  example,  day  0  in  the  model  corresponds  to  January  1,  but  day  0  in  the  migration 
pattern  corresponds  to  the  start  of  the  spring  fish  run  on  March  1),  the  offset  between  the  two 
should  be  inserted  in  days  (for  the  prior  example,  March  1  =  julian  day  60,  January  1  =  julian 
day  1,  therefore  JDO  =  (60  -  1)  =  59). 

IFRC  controls  what  fraction  of  the  population  in  the  migration  file  is  actually  run  through  the 
model.  This  card  allows  for  comparison  of  the  effects  of  differing  numbers  of  populations,  while 
maintaining  an  identical  migration  pattern.  Similarly,  FFRC  governs  the  fraction  of  the  fishing 
mortality  (set  in  a  separate  card  in  the  control  file)  to  provide  a  convenient  method  to  investigate 
the  impact  of  varying  this  parameter. 

NGRP  dictates  the  number  of  groups/classes  for  which  bioenergetics  parameters  are  to  be  pre¬ 
scribed.  IGWT  and  IGFN  are  the  individual  fish  entry  weight  and  length  for  each  class;  the 
number  of  columns  should  correspond  to  the  number  of  groups  given  in  NGRP.  EPLK,  EZOO, 
and  EDET  are  the  energy  densities  (in  J/g  of  wet  weight)  of  phytoplankton,  zooplankton,  and 
detritus,  respectively. 

The  next  set  of  cards  governs  the  bioenergetics  parameters  for  the  fish  themselves;  once  again, 
there  should  be  a  column  for  each  group  specified  in  NGRP.  These  variables  have  been  previ¬ 
ously  described  and  are  defined  in  the  variable  list,  marked  with  an  asterisk  indicating  they  are 
control  cards. 

AFNM  and  AFFM  are  the  instantaneous  annual  predation  and  fishing  mortality  rates  (fish/year), 
respectively.  The  next  sets  of  cards  (COMPOSITION  through  S  RECYC)  govern  various  para¬ 
meters  of  the  subject  fish  species;  these  cards  have  also  been  previously  described. 
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The  TRKC  card  governs  turning  the  fish  output  files  on  and  off.  Output  files  can  be  generated  at 
varying  time  intervals  (e.g.,  one  output  every  10  days,  one  output  every  365  days,  etc.);  NTRK 
governs  how  many  output  files  will  be  generated.  TRKD  governs  the  starting  day  for  each  of  the 
output  files,  and  TRKF  the  frequency  (in  julian  days)  to  output  to  each  file. 

NITKF  dictates  the  number  of  input  fish  movement  files,  which  are  listed  in  ITKFN;  these  are 
assumed  to  be  in  chronological  order  by  the  model. 

TRKFN,  MRTFN,  BINFN,  BNTFN,  and  BCONF  govern  the  file  names  for  the  various  fish  out¬ 
put  files. 

Fish  movement  is  controlled  via  fish  migration  input  files  (Appendix  D)  read  into  the  model; 
names  of  the  files  are  controlled  with  the  ITKFN  card  in  the  fish  control  file. 

The  first  line  of  the  migration  input  file  allows  the  user  to  input  information  regarding  the  run, 
and  is  not  used  by  the  model.  The  migratory  period  indicates  at  what  time  the  model  should 
move  to  the  next  fish  control  file  (if  one  exists).  The  user  should  ensure  that  the  end  time  of  the 
last  migratory  file  is  later  than  the  end  of  the  model  run  to  prevent  unpredictable  behavior  if  the 
model  attempts  to  move  past  the  end  of  the  last  input  file. 

The  rest  of  the  file  dictates  the  location  of  a  school  at  any  given  point  in  time.  At  each  time-step, 
the  number  of  schools  is  the  total  number  undergoing  movement  at  that  time-step;  any  school  not 
explicitly  listed  is  assumed  to  occupy  the  same  model  cell  as  it  did  when  last  specified.  For 
example,  school  1  enters  the  grid  in  cell  72  at  time  0.00,  and  is  assumed  to  remain  in  cell  72  until 
time-step  0.50,  when  it  moves  again.  Any  school  that  does  not  exit  before  the  end  of  the  migra¬ 
tory  period  is  assumed  to  remain  on  the  model  grid;  users  should  take  care  to  ensure  all  schools 
either  exit  the  grid  before  the  end  of  the  input  file,  or  that  their  movement  is  specified  in  subse¬ 
quent  input  files.  The  “size”  of  each  school  is  the  initial  population;  actual  population  at  any 
given  point  in  time  is  governed  by  mortality  within  the  model.  The  group  ID  corresponds  to  the 
class  number  within  the  fish  control  file  governing  the  bioenergetics  parameters. 

ADDITIONAL  INFORMATION:  For  additional  information,  contact  Dr.  P.  Soupy  Dalyander 
(601-634-4612,  Patricia.A.Dalyander@usace.army.mil )  or  Dr.  Carl  F.  Cerco  (601-634-4207, 
Carl.F. Cerco@usace.army.mil).  This  research  was  funded  through  the  System-Wide  Water 
Resources  Program  (SWWRP);  for  more  information  on  SWWRP,  please  visit 
https://swwrp.usace.army.mil,  or  contact  the  Program  Manager,  Dr.  Steven  L.  Ashby 
{Steven.L.Ashby@usace.  army.  mil). 
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APPENDIX  A:  Variables:  Variables  marked  with  an  asterisk  are  cards  in  the  fish  control  file. 


*ACTF 

*ACTNF 

*AL 

AOCR 

Bk 

*BL 

C' 

C’ 

1-  'onikM 
ConcHL 

Conphyto 

CtWP 

ctwz 

[DO] 

*DWWW 

DWWWdet 

Epred 

E0 

eff 

*EDET 

*EPLK 

*EPRD 

ZRRATE 

*EZOO 

F' 

F' 

1  c 

*FFRC 

*FTHR 

*FA 


=  activity  multiplier  during  feeding  times  (unitless) 

=  activity  multiplier  during  non-feeding  times  (unitless) 

=  empirical  fish  coefficient  for  determining  length/weight  relation  (units  may 
vary) 

=  oxygen-to-carbon  ratio  in  respiration  (gOo/gC) 

•5 

=  concentration  of  algal  group  k  (gC/m  ) 

=  empirical  fish  coefficient  for  determining  length/weight  relation  (units  may 
vary) 

=  rate  of  energy  consumption  (J/s) 

=  consumption  rate  of  element  E  (gE/s) 

=  concentration  of  detritus  (g/nr  ) 

=  concentration  of  chlorophyll  (mg/1) 

=  concentration  of  phytoplankton  (gC/m  ) 

=  ratio  of  carbon  to  wet  weight  of  phytoplankton  (gC/gWW) 

=  ratio  of  carbon  to  wet  weight  of  zooplankton  (gC/gWW) 

=  concentration  of  dissolved  oxygen  (mg/1) 

=  dry  weight  to  wet  weight  ratio  (unitless) 

=  ratio  of  dry  weight  to  wet  weight  of  detritus  (gDW/gWW) 

=  energy  density  of  the  modeled  fish  (J/g) 

=  oxycalorific  coefficient  (J/g) 

=  fish  filtration  efficiency  (unitless) 

=  energy  density  of  detritus  (J/g  of  wet  weight) 

=  energy  density  of  phytoplankton  (J/g  of  wet  weight) 

=  empirically  determined  predator  energy  density  (units  may  vary) 

=  recycling  rate  of  element  E  (g/s) 

=  energy  density  of  zooplankton  (J/g  of  wet  weight) 

=  rate  of  energy  loss  of  egestion  (J/s) 

=  egestion  rate  of  element  E  (gE/s) 

=  fish  feeding  fraction  (unitless) 

=  chlorophyll  threshold  for  fish  to  swim  at  “feeding”  speed  (g/m3) 

=  empirical  coefficient  of  egestion  (unitless) 
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*F£DW  =  fraction  of  element  E  in  ideal  fish  composition,  e.g.,  FCDW  would  denote 
fraction  of  carbon  (gE/gDW) 

*FMCF  =  recycling  of  model  constituent  MC  in  elemental  cycle  E  as  a  fraction  of 

PRRATE  from  a  living  fish  (e.g.,  fraction  of  recycled  carbon,  CRRATE, 
recycled  to  labile  dissolved  organic  carbon  is  denoted  FLDOCF) 

*FMCFM  =  recycling  of  model  constituent  MC  in  elemental  cycle  E  as  a  fraction  of  the 
total  rate  of  mass  of  element  E  returned  due  to  mortality  (see  FMCF) 

FSOE  =  concentration  of  an  organic  model  constituent,  where  F  is  the  form  (labile  or 

refractory),  S  is  the  state  (particulate  or  dissolved),  and  E  is  the  element 
(carbon,  nitrogen,  or  phosphorus);  e.g.,  RPOC  is  Refractory  Particulate 
Organic  Carbon,  LDON  is  Labile  Dissolved  Organic  Mtrogen,  etc.  (g/m  ) 

G'e  =  growth  rate  of  element  E  (gE/s) 

gap  =  fish  mouth  gap  area  (m  ) 

L  =  length  (mm) 

LZ  =  concentration  of  mesozooplankton  (gC/m  ) 

M  =  instantaneous  mortality  (fish/s) 

Nt  =  number  of  individuals  at  time  t  (#  of  fish) 

Need/,-  =  amount  of  element  E  needed  during  a  time-step  to  attain  ideal  composition 

(gE) 

*QR  =  empirical  respiration  coefficient  (units  may  vary) 

R’  =  rate  of  energy  loss  to  respiration  (J/s) 

R'c  =  respiration  rate  of  carbon  (gC/s) 

R'0  =  rate  of  oxygen  consumption  (g/s) 

*RESA  =  empirical  respiration  coefficient  (units  may  vary) 

*RESB  =  empirical  respiration  coefficient  (units  may  vary) 

*RTM  =  empirical  respiration  coefficient  (°C) 

*RTO  =  empirical  respiration  coefficient  (°C) 

S'  =  rate  of  energy  loss  to  specific  dynamic  action  (J/s) 

*SDA  =  empirical  coefficient  of  specific  dynamic  action  (unitless) 

T  =  temperature  (°C) 

u  =  fish  swimming  speed  (m/s) 

U'  =  rate  of  energy  loss  to  excretion  (J/s) 

U'c  =  excretion  rate  of  element  E  (gE/s) 

*UA  =  empirical  coefficient  of  excretion  (unitless) 
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V'  =  fish  volumetric  clearing  rate  (m3/s) 

VELNF  =  fish  swim  velocity  when  not  feeding  (m/s) 

VELF  =  fish  velocity  when  feeding  (m/s) 

WW  =  wet  weight  (g) 

AN  =  number  of  individuals  lost  within  a  time-step  (#  of  fish) 

At  =  length  of  a  time-step  (s) 
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APPENDIX  B:  Portion  of  a  water  quality  control  file,  showing  the  addition  of  the  card  to 
implement  the  fish  module. 


control  rile  ror  wqm 


TITLB  C  . TTTLBr . . 

Control  file  for  30-box  model  of  Chesapeake  Bay 
Created  Feb  23,  2005 

Consistent  with  2002  Chesapeake  Bay  model 
Seme  modifications  as  per  Patuxent  model 


These  include  flx-up  to  sed  model  CH4,  alterations  In  mineralization 
rates  and  alteration  to  algal  nitrogen  preference 

GB0M  DEFINE  NB 
4073 

NSB 

729 

NQF 

9874 

NHQF 

4530 

NSHQF 

1251 

NL 

15 

TIME  00N 

TMSTRT 

0.0 

TWEND 

1095.0 

4  DLT 

NDLT 

1 

DLT  DAY 

DLTD 

0.0 

DLTD 

DLTD 

DLTD 

DLTD 

DLTD 

DLTE> 

DLTD 

DLTD 

DLT  VAL 

DLTVAL 

864. 

DLTVAL 

DLTVAL 

DLTVAL 

DLTVAL 

DLTVAL 

DLTVAL 

DLTVAL 

DLTVAL 

DLT  MAX 

DLTMAX 

1000.0 

DLTMAX 

DLTMAX 

DLTMAX 

DLTMAX 

DLTMAX 

DLTMAX 

DLTMAX 

DLTMAX 

DLT  FTH 

DLTFTN 

0.  95 

DLTFTN 

DLTFTN 

DLTFTN 

DLTFTN 

DLTFTN 

DLTFTN 

DLTFTN 

DLTFTN 

HM  DLT 

AHMDLT 

45000 . 0 

PID3TH 

344.583 

CONTROLS 

SB  DC 
ON 

AUTOC 

ON 

VBC 

ON 

BPOC 

STLC 

ON 

ICIC 

BINARY 

ICOC 

ON 

SAVMC 

ON 

CONTROLS 

SUSFDC 

DEPFDC 

L0XC  / 

f  FISH 

ON 


ON 


ON 


2FO 

FILB . 

outputs/ 

'wqm_zfo, test_OSO506 

. . . ZP0FH . . 

BFO 

FILB . 

output s/wgm_bfo. test_050506 

. . . BF0FN .  . 

SVO 

PILE . 

outputs/ 

'wgm_s vo . te  s  t_0505  0  6 

.  . . BF0FN . . 

AHY 

FILB . 

outputs/ 

'wqm_aho . test_OSO506 

.  .  .  AHYFN .  . 

17 


ERDC  TN-SWWRP-10-1 
April  2010 

APPENDIX  C:  Example  of  a  fish  control  file. 


Control  file  for  Fish 


TITLE  C  . TITLE.  .  . 

30  Box  water  quality  model. 

Hacking  menhaden  with  30  box  sediment  transport. 

December  4,  2008 

Fish  control  file 

Extra  comment  line 

Extra  comment  line 


JULIAN 

JDO 

IFRC 

FFRC 

0 . 0 

1 . 00 

1 . 00 

GRP  PARA 

NGRP 

4 

INIT  WT 

IGWT 

IGWT 

IGWT 

IGWT 

1. 00 

85.6 

167 . 1 

284 . 8 

INIT  LTH 

IGLN 

IGLN 

IGLN 

IGLN 

48 . 0 

202.0 

252.0 

3  00.0 

ENRGY  DEN 

EPLK 

EZOO 

EDET 

602  0  . 

2  7  90  . 

900  . 

BIOE  PARA 

EFRD 

EPRD 

EPRD 

EPRD 

3  93  8  . 

8717  . 

8717. 

8717  . 

DWWW 

DWWW 

DWWW 

DWWW 

0.216 

0 . 334 

0.3  34 

0.334 

RESA 

RESA 

RESA 

RESA 

3. 

3 Ole- 3 

0.00294 

0.0027 

0.003 

RE  SB 

RES  E. 

RE  SB 

RE  SB. 

- 

■0.2246 

-0.0085 

1 

o 

o 

1 

o 

o 

RTO 

RTO 

RTO 

RTO 

33. 

33. 

33. 

33. 

RTM 

RTM 

RTM 

RTM 

36. 

36. 

36. 

36. 

QR 

QR 

QR 

QR 

2. 07 

2.5 

2.5 

2 . 5 

SBA 

SDA 

SEA 

SDA 

0.172 

0 . 172 

0.172 

0.172 

ACTF 

ACTF 

ACTF 

ACTF 

3 . 5 

3.5 

3.5 

3 . 5 

ACTNF 

ACTNF 

ACTNF 

ACTNF 

1 . 0 

1.0 

1.0 

1 . 0 
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FA 

FA 

FA 

FA 

0. 14 

0 .14 

0  .  14 

0.14 

UA 

UA 

UA 

UA 

0.1 

0.1 

0.1 

0 . 1 

MORTALITY  AFNM 

AFNM 

AFNM 

AFNM 

1.23 

0  .172 

0  .  60 

0.5  5 

AFFM 

AFFM 

AFFM 

AFFM 

0. 02 

0.22 

0  .  85 

1.37 

COMPOSITION  FCDTaT 

FNDW 

FPDW 

0.56b 

0.080 

0.024 

FSH  PAPA  GAPA 

GAPE 

SPDA 

SPDB 

SPDC 

0.02586 

1. 79767 

2.5 

-0.398 

6.378 

EFFA 

EFFB 

EFFC 

AL 

BL 

0.5 

-0.0527811  2. 

96973 

7 . IE- 6 

3 . 07 

VELF 

VELNF 

FT  HR 

FFRC 

1.  67 

0.47 

0.004 

0.75 

C  RECYC  FLDOCF 

FRDOCF 

FLPOCF 

FRPOCF 

0.25 

0.0 

0 . 50 

0.25 

FLDOCFM 

FRDOCFM 

FLPOCFM 

FRPOCFM 

0.25 

0 . 0 

0 . 5  0 

0.2  5 

N  RECYC 

FNH4F 

FUREAF 

FLDONF 

FRDONF 

FLP'ONF 

FRPONF 

0 . 5  6 

0 . 0 

0 . 24 

0 . 0 

0.1 

0 . 1 

FNH4  FM 

FUREAFM 

FLDONFM 

FRDONFM 

FLPONFM 

FRPONFM 

0 . 5  5 

0 . 0 

0.2 

0 . 0 

0.2 

0 . 05 

P  RECYC 

FP04F 

FLDOPF 

FRDOPF 

FLPOPF 

FRPOPF 

0 . 5 

0 . 4 

0 . 0 

0 . 05 

0 . 05 

FP04 FM 

FLDOPFM 

FRDOPFM 

FLPOPFM 

FRPOPFM 

0.5 

0.4 

0.0 

0.05 

0 . 05 

S  RECYC 

FPSF 

FDSF 

0.5 

0.5 

FSH  RPRT 

TRKC 

NTRK 

ON 

1 

FSH  PAY 

TRKD 

TRKD 

TRKD 

TRKD 

TRKD 

TRKD 

TRKD 

TRKD 

TRKD 

0 

FSH  FRQ 

TRKF 

TRKF 

TRKF 

TRKF 

TRKF 

TRKF 

TRKF 

TRKF 

TRKF 

1.0 

#  FILES 

NITKF 

10 
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INPUT  TRACK  FILES . ITKFN 

merged  track s.npt 
merged  track s.npt 
merged  tracks. npt 
merged  tracks. npt 
merged  track s.npt 
merged  track s.npt 
merged  track s.npt 
merged  track s.npt 
merged  tracks. npt 
merged  tracks. npt 


OUTPUT  SNAP  TRACK  FILE . TRKFN 

outputs/fsh  t tie:,  opt 

OUTPUT  MORTALITY  FILES . MRTFN 

output s / f sh  mrt . opt 

OUTPUT  MASS  BINARY  FILE . BINFN 

outputs/fsh  mass. opt 

OUTPUTS  TRACK  BINARY  FILE . BNTFN 

outputs/fsh  btrk.opt 

OUTPUTS  CONSUMPTION  BINARY  FILE . BCONF 

outputs/ consume . opt 
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APPENDIX  D:  Example  of  Migration  File. 


Atlantic  menhaden  moving  on  4,000  cell  Chesapeake  Bay  grid,  4/7/2009 
Migratory  Period 
365.25 

Time 
0.00 
#  Schools 
3 


sid 

event 

box 

size 

#  grps 

gid 

frac 

gid 

frac 

1 

Entry 

72 

6.23e+005 

1 

1 

1.00 

2 

Entry 

72 

3.81e+005 

2 

2 

0.88 

3 

0.12 

3 

Entry 

72 

7.63e+005 

1 

4 

1.00 

Time 
0.25 
#  Schools 
2 


sid 

event 

box 

size 

#  grps 

gid 

frac 

gid 

frac 

2 

Move 

71 

3.81e+005 

2 

2 

0.88 

3 

0.12 

3 

Move 

70 

7.63e+005 

1 

4 

1.00 

Time 
0.50 
#  Schools 
3 


sid 

event 

box 

size 

#  grps 

gid 

frac 

gid 

frac 

1 

Move 

71 

6.23e+005 

1 

1 

1.00 

2 

Move 

72 

3.81e+005 

2 

2 

0.88 

3 

0.12 

3 

Move 

69 

7.63e+005 

1 

4 

1.00 

Time 

1.00 

#  Schools 
2 


sid 

event 

box 

size 

#  grps 

gid 

frac 

gid 

frac 

2 

Exit 

72 

3.81e+005 

2 

2 

0.88 

3 

0.12 

3 

Move 

72 

7.63e+005 

1 

4 

1.00 

Time 
1.25 
#  Schools 
2 


sid 

event 

box 

size 

#  grps 

gid 

frac 

1 

Exit 

72 

6.23e+005 

1 

1 

1.00 

3 

Exit 

72 

7.63e+005 

1 

4 

1.00 
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